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ABSTRACT 

Detailed photoionization and radiative acceleration of self-similar magnetocentrifugal accretion disk winds 
are explored. First, a general-purpose hybrid magnetocentrifugal and radiatively-driven wind model is defined. 
Solutions are then examined to determine how radiative acceleration modifies magnetocentrifugal winds and 
how those winds can influence radiative driving in Active Galactic Nuclei (AGNs). For the models studied here, 
both radiative acceleration by bound-free ("continuum-driving") and bound-bound ("line-driving") processes 
are found to be important, although magnetic driving dominates the mass outflow rate for the Eddington ratios 
studied (L/Lndd = 0.001-0.1). The solutions show that shielding by a magnetocentrifugal wind can increase the 
efficiency of a radiatively-driven wind, and also that, within a magnetocentrifugal wind, radiative acceleration 
is sensitive to both the column in the shield, the column of the wind and the initial density at the base of the 
wind. 

Subject headings: galaxies: active — galaxies: jets — hydrodynamics — MHD — radiative transfer — quasars: 
general 
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1. INTRODUCTION 

A variety of observational signatures point to the im- 
portance of outflowing gas within many types of Active 
Galactic Nuclei (AGNs). Blueshifted absorption features 
(in Broad Absorpti on Lin e Quasars, or BALQSOs; see, 
e.g.. IWev mann et aflll991l) are seen in approximately 15% 
(Reic hardetalJl2003l) of radio-quiet quasars, with velocities 
up to 0.1c. In addition, radio-loud quasars display relativis- 
tic, collimated outflows. More recently, both UV and X-ray 
absorbing gas have been observed in ap proximately 60% of 
Seyfert 1 galaxies l lCrenshaw et alJfT999l) . Observational es- 
timates hint that the mass outflow rate is nearly equal to the 
mass inflow rate (for a review of m ass outflow in AGNs, see 
ICrenshaw. Kraemer. & Georgd2003h . 

The development of models to explain the mass outflow 
rates, geometry, and general kinematics of these winds has 
proven difficult, but progress on a number of possibilities has 
been encouraging. From very early o n, researchers exam- 
ined both radiative wind models (e.g., [Drew & Boksenberg 
Il984t IVitello & Shlosmanlll988l) and, to explain radi o jets, 
hydromagnetic models (e.g., Blandford & Pavnell 19821 here- 
after BP82). Later models of radiatively-driven winds were 
able to explain the BALQSO outflow velocities and popu- 
lation fraction ilMurrav et alj|1995l herea fter, MCGV95) as 
well as the single-peaked emission lines ( Murray & Chiang 
(19971) . The line-driven mode ls were also demonstrated 
in hydrodynamical simulations i Proga. Stone & Drewlll998t 
iProga. Stone. & Kallmanl 120001: IPerevra et alJ 12004 . which 
elucidated the density structure, geometry, and kinematics 
of the flow in two dimensions. In addition, models with 
combinations of continuum- and line-driv ing by X-rays were 
presented by Chelouche & Netzer (2001). One of the con- 
cerns with line-driving in AGNs has been the possible over- 
ionization of the gas by X-rays in the central continuum 
(leaving the wind with too few lines to intercept flux in 
the lines): the need for "shielding gas" to prevent this ove- 
rionization has been a persistent concern (e.g., M CGV95 , 
iChelouche & NetzeJ2003fllProga & Kallmanl2004l) . 
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In addition, hydromagnetic winds have also been devel- 
oped to gain insight into the observations; some of these 
models have included radiative driving. In these magneto- 
centrifugal winds, gas is loaded onto and tied to large-scale 
magnetic field lines; those field lines then fling the matter 
centrifugally away from the disk, like beads along a wire. 
Magnetocentrifugal acceleration is commonly used to ex- 
plain large-scale collimated outflows in radio galaxies (BP82; 
for reviews of m agnetocentrifugal dri ving, see ISpruitl[T99rl 
K6nig l&Pudritzl l2000: Ferre irall2003l) . In the context of the 
Broad Emission Line Region (BLR), magnetocentrifugal out- 
flows made up of clouds were first called upon to explain 
the single-peak ed broad emission li nes, outflow velocities, 
and densities (Emmering et alJ["l992l) . The "torus" (obscur- 
ing gas that plays a central role in the Unification model, 
yielding a dependence of ob served properties with inclination 
angle; see Antonucci 1993) has been explained as a dusty, 
continuous magnetic wind (Konigl & Kartie 119941 hereafter 
KK94) where radiative acceleration on dust affects the wind 
geometry. Hydromagn etic disk winds with ra diation pressure 
were also developed by de Kool & Begelmari| ill 9951) as an al- 
ternative explanation for the population fraction of BALQ- 
SOs and to understand cloud confinement within the out- 
flow. In addition, another magnetocentrifugal wind model 
was developed to explain single pe aked emission lines aris- 
ing from a distribution of clouds ( Bottorff et al] 119971) . as 
well as the dynamics of the warm absorber in NGC 5548 
(Bottorff. Korista. & Shlosman 2000). Finally, some effects 
of magnetic fields (not including magnetocentrifugal driving 
as in BP82) have also been considered in two -dimensional hy- 
drodynamic simulations (Proga 2000, 2003). 

At the present time, both radiatively driven winds and mag- 
netocentrifugally driven winds (with radiative acceleration 
added in some models) offer compelling but competing pic- 
tures of the key physics in the cores of AGNs. No clear obser- 
vational evidence yet discriminates the dominant physics of 
wind-launching in AGNs. For radiatively-driven winds, three 
main lines of evidence hint that radiative dri ving should be 
import ant in models of wind dynamics. First. ILaor & Brandl 
(2002) find a correlation between UV luminosity and the ob- 
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served outflow velocity, which agrees with a basic predic- 
tion of line-driven wi nd models dPr pga 119991) . In addition, 
the "gh ost of Ly-a" lArav et alJl!995t IChelouche & Netzerl 
2003b), as well as the realization that the radiative momen- 
tum removed from the continuum in blueshifted absorption 
lines in BALQSOs is a significant fraction of the total ra- 
diative momentum (MCGV95) are both also important clues 
that radiative acceleration is an important component to any 
self-consistent model of AGN winds. There are still, how- 
ever, concerns about how th e "shieldin g" of the win d works 
( IChelouche & Netzerl200"3btlProga & Kallmarl2004T) . In ad- 
dition, in the case of stellar disk winds, observations of two 
nova-like variables seem to s how that their winds are not dom- 
inated by radiative driving (Hartle y" st alJl2002i) . (This con- 
clusion, however, rests on the prediction that line equivalent 
widths are direct measures of mass ou tflow rate, which may 
not be the case; see Perevra et al. 2004) 

On the other hand, the leading model for collimated radio 
jets in AGNs (BP82) already calls upon large scale, dynami- 
cally important magnetic fields, as do the hvdromagnetic wind 
models (mentioned above, e.g. lKonigl & Kartielll994l|Kartiel 
119951 iBottorff et alJll997l) that have also had success in ex- 
plaining AGN observations. In addition, such a hydromag- 
netic wind would have no difficulty with overionization, and 
so might naturally serve as a "shield" for radiative accelera- 
tion further from the central source. There are, however, cur- 
rently no models which address the interplay between magne- 
tocentrifugal driving and radiative acceleration; if such wind 
models could be constrained, we may be able to observation- 
ally distinguish the physics of wind launching in the cores of 
AGNs, and gain insight into the role that outflows play both 
in accretion and in feedback of those winds in the galaxy and 
surrounding matter. 

This paper develops a detailed photoionization and dynami- 
cal model for magnetocentrifugal winds in AGNs. The model 
is designed to explore the radiative transfer within such mag- 
netocentrifugal winds, but also to help understand how radia- 
tive driving impacts the kinematic structure of such winds. 
Constructing such a model builds the foundation for later 
work to determine absorption and emission line profiles in or- 
der to compare with observations and check for the presence 
of magnetocentrifugal winds within AGNs. An overview of 
this model is presented in the model is then defined in 
detail in &\ An examination of the structure of a particular 
"fiducial" magnetocentrifugal, radiatively-accelerated wind is 
described in |4] anc l then the dependences of radiative acceler- 
ation on some initial parameters are shown in ^5] Conclusions 
and directions for future work are summarized in ^6] 

2. MODEL OVERVIEW 

Before examining the model in detail, it is instructive to 
present a summary of the basic design. The semianalytic 
model developed here includes magnetic acceleration and ra- 
diative acceleration of a continuous wind launched from an 
accretion disk. A detailed treatment of radiative transfer is in- 
cluded by using Version 96.00 of the photoionization simula- 
tion program Cloudy, last described by Ferland et al. ( 1998). 
These elements are introduced in an approximate schematic of 
the wind model shown in Figure ^ depicting a portion of the 
accretion disk and outflowing wind. In this figure, radiation 
(entering from the left side of the schematic) first encounters 
a purely magnetocentrifugally accelerated wind, which will 
be referred to as the "shield", as it absorbs radiation from 
the central continuum. The shield is introduced as a sepa- 



rate component in order to cleanly differentiate the effect of 
shielding from radiative acceleration; radiative driving of the 
shield is therefore not considered in this work. Beyond that 
shield is an optically thin, radiatively and magnetically ac- 
celerated wind streamline (which we will henceforth refer to 
as the "wind"). In this portion of the model, both magneto- 
centrifugal and radiative forces help accelerate the flow off of 
the accretion disk; the magnetic fieldlines are shown by the 
black lines bordering the outflow. The included radiative ac- 
celeration is calculated by first simulating the photoionization 
within both the shield and the wind along radial paths such as 
the thick, black lines in the figure. 

The separation of the wind into two components is, of 
course, artificial. In reality, radiative acceleration would grad- 
ually increase in importance for portions of the wind that 
are increasingly shielded. However, splitting the outflow into 
these two components allows a first-order, qualitative solution 
that can be used to gain some understanding of how magnetic 
and radiative forces might interact, and how a magnetic wind 
may be able to act as a radiative "shield" to allow more ef- 
ficient radiative acceleration. While artificial, this method of 
using two wind components has already been used success- 
fully to examine winds with magnetocentrifugal and radiative 
driving on dust (e.g., KK94, Kar tiell995l) . 

Figure|2]presents a schematic flow chart of how model cal- 
culations proceed. The wind starts as a self-similar magneto- 
hydrodynamic model that yield s the pure magnetocentrifugal 
wind solution (covered in 33. 11 1. Next, simulatio ns of the pho- 
toionization balance of that wind streamline ( §3.2> are run, 
and the resultant ionization balance and transmitted contin- 
uum are us ed to calculate the radiative acceleration behind 
the shield ( 33. 3> . Next, the radiative acceleration is input (as 
a function of polar angle, 9) back into the self-similar magne- 
tohydrodynamic model, modifying the structure of the wind 
streamline, while leaving the shield unaffected. This process 
is then repeated, simulating the photoionization of that mod- 
ified wind and recalculating the radiative acceleration terms. 
We typically iterate five to eight times to converge to a final 
equilibrium solution. 

With the basic model now summarized, it may be instruc- 
tive to com pare and con trast it to the recent wind model 
examined in Proga (2003), where the combination of mag- 
netic and ra diative forces in disk winds is also investigated. 
Proga (2003) concentrates on numerical simulations of time- 
dependent winds with line driving and magnetic forces. These 
numerical simulations allow large-scale models of outflows 
that are valuable in understanding global wind structures 
in many different astrophysical contexts. In contrast, the 
radiatively-driven components of this model are more local- 
ized, since radiative acceleration is applied to the shielded 
streamline only. This model setup has been chosen for its 
flexibility in radiative acceleration modeling; using Cloudy 
for such radiative acceleration models enables computations 
not only of radiative line driving but also of continuum driv- 
ing, and allows us the freedom to include dust and easily vary 
the incident spectru m, for exam ple. In addition, the magnetic 
winds produced in Proga (2003) are not magnetocentrifugal 
outflows (as in BP82), as are the semianalytic winds presented 
here. Further, these new models are steady-state, not time- 
dependent. Finally, semi-analytic, steady-state models allow 
an exploration of general behaviors through many parameter 
variations; large-scale numerical simulations can usually vary 
only a few parameters. In summary, these models cover dif- 
ferent facets of the disk-wind problem, yielding different per- 
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FIG. 1. — Schematic of the basic geometry and major components of the radiatively accelerated magnetocentrifugal wind model. The wind itself is split 
into two components: a pure magnetocentrifugal wind that acts as a shield (the wider component on the left) and an optically thin streamline with combined 
magnetocentrifugal and radiative acceleration. The heavy black lines indicate a few of the radial zones where Cloudy simulates the photoionization balance of 
the wind (typically, ~ 40 such radial Cloudy simulations are run, spaced logarithmically in polar angle, 9). 
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FIG. 2. — The iterative scheme in the wind model, showing a single iteration of the model, moving counterclockwise around the figure. The model starts by 
solving the self-similar magnetohydrodynamic (MHD) equations for the structure of a centrifugally driven outflow (upper left). Next, photoionization simulations 
determine the continuum incident on the radiatively accelerated streamline as well as the ionization state of the gas in the streamline, and then those results are 
utilized to calculate the radiative acceleration of that flow. The radiative acceleration is then applied (as a function of the polar angle 9) to the second component 
of the self-similar magnetocentrifugal wind. A new wind structure for that second section of the outflow is then derived from the interaction of those forces 
(upper right). 



spectives on a complicated system. 

3. THE TWO-PHASE HYDROMAGNETIC AND RADIATIVE WIND 
MODEL 

In this section, we describe in detail the model's com- 
ponents, and derive key equatio ns. T he magnetocentrifugal 
wind model is introduced first ( 33. 1> . followed by the pho- 
toionization simulations ( 33. 2> . radiative acceleration calcu- 
lation ( 33.3> . an d fina lly a discussion of the wind model equa- 
tion of motion ( 33.4V 

3.1. Magnetocentrifugal Self-Similar Wind Solution 



To derive the equations governing the continuous magneto- 
centrifugal wind, we start with the equations of a stationary, 
axisymmetric magnetohydrodynamic (MHD) flow in cylin- 
drical coordinates, make the assumption of self-similarity 
in the spherical radial coordinate, and utilize the continuity 
equation, conservation of angular momentum along the flow, 
and both the radial and vertical momentum equations, very 
much as in BP82 and KK94 (see Appendix lAV Thermal ef- 
fects are neglected in the wind, therefore effectively assuming 
that the wind start s out supersonic (this assumption is checked 
later on, see ^3.51 . In deriving the wind equations, we use the 
same simplifications as BP82, except for the added complica- 
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tion that energy is not conserved in the radiatively-accelerated 
system due to the constant input of radiative energy into the 
wind. In the original formulation of BP82, conservation of 
energy supplied an additional constraint which allowed a sim- 
plification of the equations of motion to two first-order differ- 
ential equations. Because energy is not conserved in this flow, 
the equivalent of three first-order differential equations must 
be integrated, solving for three parameters simultaneously in- 
stead of two as in BP82. The detailed setup and derivation of 
this set of equations of motion are given in AppendixlAl 

The integration of the momentum equations starts by spec- 
ifying the following initial parameters: the mass loading of 
the wind (the ratio of mass flux to magnetic flux in the mag- 
netocentrifugal wind, k oc 47r f Vp , where p is the mass den- 
sity of the wind, v p is the poloidal velocity of the wind 
(Vp = (v^ + Vj) 1 / 2 ), and B p is the poloidal magnetic field 
strength), the specific angular momentum of gas and field in 
the wind, and the power-law exponent (b) that describes the 
change in density with spherical radius: p oc R~ . Also in- 
put, as parameters, are the mass of the central black hole, M., 
the wind's launch radius on the disk, ro, and the density at 
the base of the wind at the launch radius, no- The program 
employs a "shooting" algorithm (using the SLATEC routine 
DNSQ; Powell 1970) to integrate from the singular point (the 
Alfven point) to the disk, solving for the height of the singu- 
lar point above the disk {\a) and the slope of the streamline at 
both the disk and the Alfven point (£' and £' A ) by matching the 
integration results to boundary conditions on the disk. After 
solving for the position of the Alfven point, the equations of 
motion are integrated from the disk to a user-specified height 
beyond the Alfven point; along the streamline, the run of ve- 
locity, density, and magnetic field are calculated. 

This code has been tested (without radiative acceleration) 
against the solution given in BP82 and have duplicated their 
results to within 8%. This is close to the previously reported 
4% va riance in the recalculation of BP82 reported in ISafierl 
( 119931) . The difference in these new results comes not only 
from using higher precision calculations compared to BP82 
(we use the same precision as Safi erll993l) . but in addition, a 
more complex set of equations is being solved. 

3.2. Photoionization Simulations of the Wind 

Next, Version 96 .00 of the photoionization code Cloudy 
dFerland et all 19981) is used to simulate the absorption of the 
magnetocentrifugal shield and wind as well as the ionization 
state at the wind streamline. Photoionization simulations of 
the shield and wind are used to calculate the radiative accel- 
eration and allow considerable flexibility in gas parameters 
(such as gas abundances, dust, central continuum, etc). This 
flexibility is gained at the cost of simulating the photoioniza- 
tion state of the wind as if it were a static medium, as Cloudy 
assumes; this is only true of the recombination timescale for 
the gas is much shorter than the transit timescale of gas in 
the region simulated (T reC omb < ifl ow ). We have however, veri- 
fied, a posteriori, that T recom b < Tfl ow for all of the radiatively 
accelerated ions at the base of the wind where radiative ac- 
celeration i s important, using the recombinati on rate approxi- 
mation s in lArnaud & Raymond! JT9 92) and Verner & Ferland 
(1996) to calculate the recombination times. Even at high lat- 
itude, Tjecomb < Tfl ow for all of the significant ions. At the high- 
est latitudes, highly ionized O, which has the longest recom- 
bination time among the significantly radiatively accelerated 
ions, has a recombination timescale that is still a factor of two 



less than the transit time. Meanwhile, highly ionized Fe ions, 
which dominate the low level of radiative acceleration at high 
latitudes, have recombination timescales an order of magni- 
tude less than the transit timescale. Using Cloudy is there- 
fore a reasonable approximation for the radiative equilibrium 
within our winds (especially since, with ou r code, adiabatic 
and advective effects are added, see 33.2. 1> . As Cloudy en- 
ables a flexible, self-consistent calculation of radiative accel- 
eration, we accept this approximation to enable these calcula- 
tions. 

Of basic importance to the photoionization simulations 
is the illuminating spectral energy distribution (SED). For 
the purposes o f these simulations, an SED adapted from 
Risaliti & Elvis (2004) is used (an example of the SED is 
shown in Figure [6ji. This SED is input into Cloudy via 
the generic "AGN" continuum with Thiackhoriv = 1.5 x 10 5 K, 
a ox = -1 .43 dElvis. Risaliti. & Zamo rani 2002), a uv = -0.44, 
and ax = 0.9 (a ox defines a single power-law that would de- 
scribe the continuum between 2500 A and 2 keV, auv is the 
slope of the low-energy component of the Big Blue Bump, 
and ax is the X-ray power-law exponent; our value of a ox 
is taken from the m iddle of the range 0.8 to 1.0 given in 
Risaliti &Elvis1 ( 12001 1. 

Spectral signatures and radiative acceleration also depend, 
of course, on the column density in the wind. Observations 
yield only rough constraints for this, so these columns as left 
as free parameters; the effect of varying these columns will be 
investigated in this paper. The columns throughout the shield 
and wind are set by the columns at the base of the shield and 
wind, denoted A^.o for the hydrogen column at the base of 
the wind. As the wind rises above the disk and accelerates, 
that column density (A^h) drops as a function of height due 
to mass conservation (an example of this is shown later in 
Fig Q- Investigating the shielding ability of such a dynamic 
shielding column is of central interest to this paper, and will 
be addressed in 35.11 

3.2.1. Photoionization of the Continuous Wind 

As depicted in Figure^ Cloudy simulates the photoioniza- 
tion of the wind along radial sight lines through the shield 
and through the wind, ending at the site of radiative accel- 
eration on the wind streamline. The continuum incident on 
the wind streamline is also calculated. The photoionization 
state and continuum at the end of the Cloudy calculations are 
recorded, and then used to compute the radiative acceleration 
of that gas. Finally, the acceleration is tabulated and applied 
as a function of 6 along the wind streamline by inputing the 
angle-dependent radiative acceleration into the gravitational 
term (this is covered in more detail in Appendix A, see eqn. 
[Alt . 

Whereas Cloudy is designed to simulate the photoioniza- 
tion balance of gases as in the shield and wind, it cannot eas- 
ily incorporate adiabatic and advection effects: Cloudy has 
no knowledge of the particular velocity profile of the wind 
in the overarching model, or the temperature difference be- 
tween successive photoionization models as the wind climbs 
above the disk. Therefore, this model calculates both advec- 
tive heating and adiabatic cooling in the wind, and adds those 
terms manually into Cloudy's simulations. Both terms largely 
cancel in the wind, and have only a negligible effect on out- 
flow dynamics, but they are included in all of the models for 
completeness. 

3.3. Radiative Acceleration Calculations 
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The model then incorporates the above-mentioned results 
for the ionization structure and radiation field to calculate the 
radiative forces felt by the wind. There are two different kinds 
of radiative acceleration to consider: continuum acceleration 
(including radiative acceleration on dust) and line accelera- 
tion. It is convenient to express the radiative acceleration in 
terms of r(6>), 



T(0): 



^radiative(^) 



(i) 



where a ra diative is the acceleration due to radiation, and g is the 
local gravitational acceleration. 

3.3.1. Line and Continuum Acceleration 

In general, for continuum and line acceleration, the radia- 
tive acceleration is given by 



r= 



pc 



(M com +M lines ) 



CM. 



(2) 



where F is the local flux (the flux transmitted through both 
the shield and wind column), n e is the electron density, p is 
the gas density, c is the speed of light, G is the gravitational 
constant, M, is the mass of the central black hole, r and z are 
cylindrical coordinates centered on the black hole, and M cont 
& Mines are the "force multipliers" that relate how much the 
radiative forces on the gas (on the line and continuum opacity, 
respectively) exceed the radiative forces on electrons alone. 
They are given below in terms of the continuum opacity, \v, 
and the line opacity, \h for the continuum and lines, respec- 
tively: 



M com = — — - \ XvF v dv, 
n e a T F J 

1^ l-e"'"' 
Mi ines = -> F,Av, . 



with 



m 



ajn e 



dvR I dR ' 



(3) 
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where v is the photon frequency, F/ is the local (transmit- 
ted) flux in the line at the frequency of line number /, v tn 
is the sound speed in the gas, Ais/ = vv^/c is the thermal 
line width, and r\i compares the opacity of the line (for a 
given ionization state of the gas) to the electron opacity, rep- 
resenting all of the atomic physics in the radiative acceler- 
ation calculation. The last remaining variable, f, is often 
called the "effective electron optical depth" and encodes the 
dynamical information of the wind in the radiative accel- 
eration calculation. This dynamical information is impor- 
tant because in an accelerating medium, one must also ac- 
count for the Doppler shift of the atomic line absorption en- 
ergy in the accelerating gas relative to the emitted line pho- 
ton's energy: beyond the Sobolev length, v t h/(dvR/dR), in- 
cluded in t, a line photon will be Doppler-shifted out of 
the thermal width of the absorption line and can escape 
the gas (see LSobolevI 119581 (C astor. Abbot t. & Kleinl 119751 
Mihalas & Weibel-Mihalasl 19991) . 

The above-mentioned force multipli ers are cal culated us- 
ing the resonance line data of IVerner et all (|1996) with solar 
abundances. Since M cont , the continuum multiplier, depends 
only on the ionization state, it is tabulated solely as a function 
of height in the wind. In contrast to M cont , the line multiplier 



(Mine) is tabulated for a range of values of the parameter t. 
Later, when calculating the equation of motion for the wind, 
the local velocity gradient is used to compute the actual value 
of t, which is then used to linearly interpolate the table of Mi; ne 
and then evaluate the radiative acceleration. 

The force mult iplier computation has been tested against 
lArav et al . (1994), who also calculated radiative acceleration 
from photoionization simulations. Figure [5] compares these 
new calculations results against their fits (noting that there is 
a typo in their eq. [2.9]; Z.-Y. Li, personal communication), 
where the force multipliers as a function of the ionization pa- 
rameter U is presented (U is the ratio of hydrogen-ionizing 
photon density to hydrogen number density n, given by U = 
Q/4TrnR 2 c, where Q is the number of incident hydrogen- 
ionizing photons per second, and R is the distance from the 
continuum source). O verall, good agreement is found, espe- 
cially considering that Arav et al. ( 1994) point out that their 
fit deviates from their calculations at low values of U. The 
incr ease in the newly -calculated continuum force multiplier 
over lArav et all (119941) is most likely due to the different con- 
tinuum opacity database included in Clou dy 96 compared to 
the code (MAPPINGS) that was used in lArav et ail ( 11994 . 
The multiplier values and trends with ionization parameter are 
still clearly very similar, however. 

3.3.2. Non-Sobolev Effects 

Simply using the Sobolev length, v t h/{dvn/dR), in Equation 
|4]can be misleading. In early simulations, we found that this 
Sobolev length could, in regions where the gas is slowly ac- 
celerating, be much larger than the physical size of the shield 
and wind combined. This is clearly not physical, so a sim- 
ple non-Sobolev method was employed to calculate the size 
of the combined shield and wind column. 

First, for the wind, the length of the absorbing column is 
simply limited to the minimum of the wind's Sobolev length 
and its true (physical) length. 

Second, for the shield, the length of the absorbing column 
is given by the minimum of the shield's Sobolev length 1 , its 
true physical length, and the average length of the column 
that absorbs photons at the wavelength of the dominant ac- 
celerating atomic lines. This last length-scale requires more 
explanation. To calculate this length, the top 20 line force 
multipliers (for individual atomic lines) are found for each 
polar angle 6, and for each of those lines, the column of each 
particular ion in the shield is read from the Cloudy simula- 
tions. An average shielding column for all of the high-opacity 
lines is then calculated. To test this approximation, we have 
used lists of the top 10, 20, and 50 transitions in the wind 
and used them to calculate the limiting column in the shield. 
Changing the number of lines included does not significantly 
change the final wind solutions, especially since this effect is 
only important at the extreme base of the wind. Without con- 
sidering all of these constraints on the size of the absorbing 
column, the Sobolev length can significantly overestimate the 
optical depth in the shield (even overestimating the physical 
length of the columns, for small accelerations), which results 
in the line driving acceleration dropping below the continuum 
acceleration. Other non-Sobolev effects (such as line blan- 
keting within the shielding gas or the wind) are not consid- 

1 The Sobolev length calculation for the shield does take into account the 
offset in velocity between the wind and shield, which is important since the 
wind is radiatively accelerated in addition to the magnetic acceleration that 
the wind and shield share. 
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FIG. 3. — Comparison between the force multiplier calculation in this code (solid lines) and the fits of Arav et al. 1 1994, dashed lines). The plot in the upper 
left is the comparison of the continuum force multiplier, as a function of the ionization parameter, U. The next three panels show a comparison of the line force 
multiplier (which is also a function of t, besides U) as a function of U for three different values of f, given at the top of each panel. 

ered in this model; for a co nsideration of these effects, see 
IChelouche & NetzeiU2003al) . 

With a length-scale found from the sum of the wind and 
shield lengths found above, that length-scale is then substi- 
tuted for the Sobolev length in the calculation of Mi; nes in 
Eq. |4] Including this physical length of the wind and shield 
will introduce dependences on the sizes of the wind and shield 
columns, which will be examined in |5] 



As we are interested in steady-state winds, ^ in Eq.[6]is set 
to zero. As already mentioned, gravitational, radiation, and 
Lorentz forces are included, while thermal terms are not. This 
yields the expression below, with the magnetic force term is 
split into pressure and tension components. 



p(vV)v=-[l-T(6)] 



GMp 



1 



1 



R Vfi z + — (B-V)B (7) 

87T 4tt 



3.4. Integrating the Euler Equation for the Wind 

Given the results from the ionization and radiative accel- 
eration calculations, the next step is to solve for the effect 
of radiative acceleration on the wind, taking the magnetocen- 
trifugal wind model already computed and augmenting its ac- 
celeration with radiative forces. To do this, the full equation of 
motion (the Euler equation) is integrated along the streamline 
of the self-similar wind solution, recording T(9), the radiative 
acceleration. T(8) is then input into the self-similar model in 
the subsequent iteration (see Appendix A; specifically, eqn. 
lA3t . 

In its simplest form, Euler's Equation is given by 

p(5 +(v - V)v )=E F <'< ( 6 ) 

where the F, represent the various forces included. 



For these calculations, it is more intuitive to integrate the 
equation of motion along the flow already given by the mag- 
netocentrifugal wind solution. Therefore, taking the dot prod- 
uct of Euler's equation with s, which is defined as the direction 
along the flow, and expanding and simplifying the left-hand 
side of the equation, one finds: 



(v • Vv) • s = v p — - - — sin B F , 
as r 



where 



: tan 



dr\ 
~dz) 



In the same way, if is defined via 



= tan 



(8) 



(9) 



(10) 
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the gravitational term can be written as 

-[l-T(ff)]^-R-S=-[l-T(6)]-^2CO & (e-e P ). (11) 
R 2 - r z +z 

Next, we take the dot product of s with the magnetic terms to 
find 



1 



VB 2 + 



1 

47T/9 



(B • V)B 



Btf, djrB^) 
47r pr ds 



(12) 



Combining all of those terms, the full Euler Equation is ob- 
tained: 



dv p 



vi GM 
r {_r l +z l ) 



cos(6-6 F ) 



djrB^) 
4npr ds 



(13) 



This equation is still dependent on v^, however, which can 
be eliminated by appealing to the induction equati on and 
I i for such axisymmetric systems (see, e.g., Konigl & Pudritz 
2000), to yield a relation between v p and v^: 



v p B^ 

B n 



+ Qr. 



(14) 



Substituting this expression into the Euler Equation yields: 



dv p 



v p B<{> 



Sinf/r 



-[i-rm 



GM 

(r 2 +z 2 ) 



COS(0-0jO- 



B<p djrB^) 
Aitpr ds 



(15) 



To evaluate the effective optical depth f , an expression is re- 
quired for dvfi/dR, the spherical radial gradient of the spher- 
ical radial velocity. Since the code integrates quantities only 
parallel to the flow, approximations to the perpendicular ve- 
locity gradients must be used. Assuming that the derivatives 
of 9 and 9 F with distance along the streamline are small (ver- 
ified a posteriori to be true): 



dv R 
dR 



= RVv R , 

= R-\7[v p cos(9-9 F )], 

«cos(6>-6V) | — sin (9 + 
V dr 



dvp 
dz 



■ COS( 



(16) 
(17) 

(18) 



( cos(6» - 9 F ) ^ (sin 9 F sin 9 + cos 9 F cos 9) , (19) 
ds 



t dv n 
= cos 2 (9-9 F )—^-. 

ds 



(20) 



This integration procedure has been tested with radiative 
acceleration turned off, where it reproduces the original self- 
similar velocity profile to within one part in 10 5 . With the ra- 
diative acceleration turned on, the entire code has repeatedly 
converged within approximately eight iterations to an equilib- 
rium magnetic wind structure (see Fig. II 1> . These tests show 
that consistent solutions are found and that radiation pressure 
does indeed affect the shielded component of the wind. 

3.4.1. Critical Points 

As with any steady-state wind dynamics problem, one 
must search for and consistently pass all critical points (e.g., 
IVlahakis et alJ 12000). Critical points mark the location in 



the wind where the flow speed is equal to the speed of in- 
formation propagation in the wind, and mark locations in 
the solutions where solutions branches, or roots, meet; in 
the case of radiatively accelerated winds, the location of the 
critical point is set by the information propagation sp eed of 
radiative-acoustic wave, or Abbott spee d (see, e.g., [Abbott 
119801 iMihalas & Weibel-Mihalasl l!999). In integrating the 
equation of motion for the radiatively-accelerated wind, this 
code searches for critical points by looking for multiple roots 
in the solution to the equation of motion. However, no critical 
points due to radiative acceleration are present in any of the 
wind solutions we have found (the magnetocentrifugal wind 
does, however, always pass through its own Alfven critical 
point). 

To check this result, we have duplicated the work of 
Feldmeier & Shlosman ( 1999, hereafter, FS99), verifying that 
for simple wind geometries and without magnetocentrifugal 
acceleration, the integration code does indeed encounter a ra- 
diative critical point as predicted and found by FS99. In par- 
ticular, for the field geometries and forces used in FS99, we 
have found identical solutions to both their analytical and nu- 
merical calculations. 

We then gradually add, to the FS99 model, new components 
that are present in our new calculations. When the centrifu- 
gal acceleration and the enforced corotation near the base of 
the wind are introduced into the framework of FS99, the cen- 
trifugal acceleration overwhelms the acceleration of the sec- 
ond root that was present in FS99. Therefore, only a single 
root is found, and with only a single root, a critical point can- 
not be present in our solutions. Another way to check this 
is to examine the limit of winds launched at very large an- 
gles to the accretion disk. Indeed, for those large angles, the 
radiative acceleration begins to dominate the centrifugal ac- 
celeration, and the critical point reappears. Therefore, for the 
geometry of these magnetocentrifugal winds, where the angle 
of the outflow to the accretion disk surface is less than 60°, 
no radiative critical point is expected within these solutions: 
a radiative critical point will not be present when centrifugal 
acceleration is dominant. 

3.5. Model Assumptions and Limitations 

This model includes simplifying assumptions about the out- 
flow in order to make these calculations possible. In this sub- 
section, the assumptions and limitations of the model are sum- 
marized, as are the reasons for allowing those assumptions. 

First and foremost among these assumptions is self- 
similarity. While enabling a relatively quick and flexible 
model that can be used to survey a wide variety of outflows, 
this assumption does impose constraints on the dynamics of 
the model. However, the assumption of self-similarity is es- 
sential to the magnetocentrifugal wind solution, as it simpli- 
fies the complicated MHD equations. Simultaneously, the ra- 
dial self-similarity accommodates, very easily, the radial ge- 
ometry for the photoionization simulations. 

The above-mentioned Cloudy photoionization simulations 
assume a static med ium, which is an approximation as well. 
As we show in 33.21 however, this approximation is valid for 
the these wind models. 

Since the accretion disk is a boundary condition in these 
models, these winds are assumed to be loaded with matter 
from the disk. Accurate models of the accretion disk struc- 
ture are beyond the scope of this paper, so it is assumed that 
the full mass outflow rate of the the wind is indeed input onto 
the magnetic fieldlines at the accretion disk surface. In ad- 
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dition, the matter that is loaded onto those fieldlines is as- 
sumed to flow supersonically, i.e., the gas has already passed 
the sonic point in the flow. This is done to simplify the mag- 
netocentrifugal wind equations and retain the basic model as 
outlined in BP82; as such, the same asymptotic expansions 
near the disk (from BP82) are utilized here. This treatment 
of the wind has been checked in several different ways. First, 
the magnetic pressure in the wind is indeed greater than the 
thermal pressure throughout the entire wind. Also, the wind's 
final velocities are much greater than the sound speed at the 
base of the wind, showing that thermal effects are negligible 
in determining the final wind velocities. Finally, the lowest 
speeds found in the wind model are of order 20% of the sound 
speed, and such low Mach numbers are found only very near 
the disk, at the base of the wind. Given the above evidence of 
the dominance of magnetic fields and radiative acceleration, 
the "cold-wind" approximation is valid for these calculations. 

Also, in the calculation of the radiative acceleration of the 
wind, an approximation to the velocity gradient along spheri- 
cal rays is required. This approximation is necessary because 
the Euler integration for the wind yields velocity gradients 
only along the flow (the poloidal velocity gradients), and thus 
the o ther components must be approximated geometrically 
(see El. 

The approximations and limitations outlined above do con- 
strain the use of this wind model, but in making these com- 
promises, a versatile tool can be developed to study the chosen 
geometries and forces. 

4. RADIATION TRANSFER WITHIN A FIDUCIAL 
MAGNETOCENTRIFUGAL WIND 

For definitiveness, this model is first employed to examine 
the radiative transfer within one fiducial irradiated magneto- 
centrifugal wind. For the purpose of this paper, 'fiducial' is 
defined to indicate the parameters listed in Tabled These pa- 
rameters are not meant to represent a proposed model for any 
one particular AGN, but to define a starting point from which 
to examine the structure of these outflows as well as the de- 
pendences of the outflows on the model parameters. For in- 
stance, the shielding column of A^H,shieid,o = 

10 23 cm (where 

A^n.shieid.o represents the value of Mi shield at the base of the 
wind, i.e. just above the accretion disk surface) is chosen be- 
cause it displays an amount of radiative acceleration between 
the extremes of the smaller and larger columns that will be 
tested. Similarly, the radiative wind column that is defined is 
again between the extremes of nearly optically thick A^H,rad,o = 
10 23 cm" 2 and very optically thin Nn,md,o = 10 19 cm" 2 . Study- 
ing such a model first will help bring into focus important 
issues concerning the interplay of dynamics and photoioniza- 
tion, and represents a foundation from which one can explore 
the parameter dependencies of the model. 

This fiducial model was therefore run with the parameters 
given in Table [2 and after eight iterations, converged to the 
final wind structure. The results for the fiducial model are 
shown in Figures 0] though In these figures, an overview 
of the equilibrium state of this model is presented. The results 
displayed in these plots are discussed in detail below. 

First, Figure @] shows the height of the poloidal streamline 
as a function of radius in units of the launching radius. In ad- 
dition, to illustrate the small difference in geometry between 
the final and initial wind models, Figure|5]gives the fractional 
change in height as a function of radius. Both of these figures 
show that the wind still maintains a collimated state, achiev- 
ing a height of z/ro ~ 100 (where ro is the launching radius) 



at a cylindrical radius of only r/rn ~ 30. So, in this fiducial 
model, despite the input from radiative acceleration, the wind 
maintains this streamline with only small changes in the struc- 
ture of the wind throughout all iterations (see Fig. [3}. Thus, 
for the case of the fiducial model with L/Lem = 0.01, this 
added acceleration does not significantly affect the structure 
of the magnetocentrifugal outflow. These models do show 
changes in the velocity structure of the wind near the disk 
surface (as will be shown in Fig. II It . but the poloidal wind 
structure does not change significantly: on the scale of Fig- 
ure|4] the streamlines of the initial, purely magnetocentrifugal 
streamline would lie on top of the streamline shown. 




r [r ] 

FIG. 4. — Poloidal wind streamlines in the fiducial model. Both the cylin- 
drical radial coordinate, r, and the height, z, are given in terms of the launch 
radius, ro. Note that the wind is still somewhat collimated, since the radiative 
acceleration input is fairly low for L/Ledd = 0.01. 
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FIG. 5. — Fractional difference in streamline height as a function of radius 
between the final iteration and first iteration of the fiducial model. The small 
different in height emphasizes that the structure of the magnetocentrifugal 
wind is not significantly modified for the fiducial model's low L/LEdd = 0.01. 
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TABLE 1 

Parameters adopted for the 'fiducial' model in this study. 



PciTum c tcr 


Fiducial Value 


Parameter Description 


no 


10 9 cm" 3 


initial density of the wind at the launch radius 


A'H.shield.O 


10 23 cm -2 


gas shielding column at the base of the wind 


iVH,rad,0 


10 21 cm" 2 


gas column, behind the shield, that is radiatively accelerated 


M m 


10 8 M© 


mass of the central black hole 


Incident Spectrum 


Risaliti & Elvis (2004) 


Spectrum for the central continuum 


^continuum 


0.01 L Edd 


luminosity of the central continuum 


'"0 


3 X 10 16 cm 


launch radius of the wind 


K 


0.03 


dimensionless ratio of mass flux to magnetic flux in the wind 


X 


30.0 


normalized total specific angular momentum of the wind 


b 


1.5 


power-law describing variation of density with spherical radius in the wind: n <x R~ b 


Dust in Wind 


No 


presence of dust in the wind 



At logarithmically-spaced co-latitudinal angles along the 
streamline, Cloudy photoionization simulations are run to de- 
termine the photoionization state of the gas, as well as the 
radiative transfer through the shield and wind. Changes in 
the continuum transmitted through the shield are shown in 
Figure |6j the various plots show the simulated continuum at 
various heights in the shield corresponding to the indicated 
columns. As the shielding column decreases as a function 
of height above the disk, the shield transmits progressively 
more and more of the ionizing radiation. This plot also dis- 
plays how rapidly the column drops as a function of height 
above the disk: A^H.shieid.o = 10 23 cm" 2 occurs at 9 = 89.9°, 
Akshieid = 10 225 cm" 2 at 9 = 89.8°, Nb^uu = 10 22 cm" 2 at 
9 = 89.4°, and A^hieid drops to 10 21 cm" 2 at 9 = 85.1°. 

The flux transmitted through the shield then illuminates the 
radiatively accelerated wind; results from the photoionization 
simulations for the wind are presented in Figure where the 
streamline, velocity, density, ionization parameter and tem- 
perature in the wind are plotted. In Figure^, the height of a 
wind streamline as a function of distance along the streamline 
(labeled s, given in units of the initial radius, ro) is shown; this 
plot simply recasts the structure of the flowline shown in Fig- 
ure 0]in terms of s for comparison with the remaining plots. 

In Figure 03, velocities along the streamline are plotted, 
showing not only the rapid acceleration in the wind, but also 
comparing the components of the wind's velocities. All ve- 
locities are given in units of the Keplerian velocity at the base 
of the wind, v^o (Vk,o = 6.65 x 10 3 km s" 1 for the parameter 
values in Table[Q. This plot shows that the vertical velocity is 
dominant in these winds at large distances (again showing the 
wind is somewhat collimated), with the radial and azimuthal 
velocities becoming approximately equal far from the launch- 
ing radius of the outflow. Near the very base of the disk, the 
radial velocity quickly dominates both the azimuthal and ver- 
tical velocities (qualitatively similar to the velocity structure 
calculated for radiatively-dominated flows as in MCGV95). 
Most importantly, though, we note the extraordinarily rapid 
acceleration of the gas from the disk; such acceleration is a 
hallmark of both magnetocentrifugal as well as radiatively- 
dominated winds, which usually accelerate to their terminal 
velocities in a distance on the order of their launching radius. 

Due to mass conservation, the extremely rapid accelera- 
tion of the wind causes a sharp drop in both the number den- 
sity and the column density with height above the disk: both 
the number density and column density immediately drop by 
three orders of magnitude as the wind rises above the disk. 



This is displayed in Figure 0;. This overall drop in density 
is extremely important for the ionization state of the wind, 
not only for the observational ramifications (i.e., what ions 
are present in various parts of the wind) but for the accel- 
eration of the wind as well, as will be shown very shortly. 
Figures |8]and|9]show in more detail how the radiative accel- 
eration leads to a substantial change (by approximately a fac- 
tor of two) in both velocity and density with height near the 
disk surface. The changes in velocity for the pure magneto- 
centrifugal wind as compared to the magnetocentrifugal and 
radiatively-accelerated wind is shown in Figure [8] The dif- 
ference in the density profile for the pure magnetocentrifugal 
wind as compared to the magnetocentrifugal and radiatively- 
accelerated wind is shown in Figure[9] (Note that all densities 
are normalized to their value at the base of the wind.) 

The corresponding ionization state of the wind and temper- 
ature are shown in Figure 01. Most striking is the dramatic 
rise in the ionization parameter as the wind rises above the 
disk, which is simply due to the drop in density and in shield- 
ing already mentioned. The ionization parameter is of prime 
interest, as the radiative acceleration in resonant lines is de- 
pendent on the number of atomic lines in the gas; the rapid 
ionization of the gas prompts questions about how efficient 
line-driving will be within this magnetocentrifugal wind. In 
addition, can a wind with such a dramatic drop in column 
density form an effective "shield"? And for these models, 
how does the radiative acceleration then compare to magnetic 
acceleration? 

These questions are addressed in Figure [H)| From the 
Cloudy simulations summarized in Figure both the bound- 
free and bound-bound radiative acceleration are calculated; 
the resultant acceleration (compared to the local gravity) for 
the fiducial model is shown in Figure ^3 This model shows 
that both line-driving and continuum-driving have important 
roles to play, with the line-driving dominating the contin- 
uum driving in the high-density, low-ionization part of the 
wind, and continuum driving dominating line-driving at larger 
distances, when the density is much lower. It is important 
to note that for the parameters of the fiducial model with 
L I ^Edd = 0.01, magnetic acceleration is still much greater than 
either line or continuum acceleration. Line-driving is greater 
than continuum-driving in only part of the outflow (although 
this can chan ge wi th the density at the base of the wind, as will 
be shown in 35. 2L as the ionization parameter is low enough 
only at the base of the wind for significant numbers of atomic 
lines to exist. In addition, radiative driving is not immedi- 
ately important at the extreme base of the outflow, because of 
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FIG. 6. — Both the incident continuum on the shield and the transmitted continuum through the shield is displayed. Note the rapid decline in the shield's 
column density as a function of polar angle: A'H.shield.o = 10 23 cm" 2 lies at the polar angle 9 = 89.9°, Nh, shield = 10 22,5 cm -2 at 9 = 89.8°, JVh, shield = 10 22 cm -2 at 
9 = 89.4°; the radiation encounters Wh. shield 10 21 cm" 2 at only 9 = 85.1°. Also indicated for each column is the distance along the flowline is) where that column 
occurs, given both in cm and in units of the launching radius, rn. 



the low fluxes that penetrate the columns there. Meanwhile, 
the acceleration due to continuum-driving stays close to the 
Eddington ratio, at ~ 0.01. This value is reasonable, as most 
of the continuum acceleration comes from electron scattering, 
so that r ~ 0.01 would be expected. Minor increases above 
that value very near the disk surface are due to bound-free 
transitions in the portion of the wind closer to the disk, where 
Tcontinuum rises to ~ 0.02. 

As has already been shown, at this low Eddington ratio the 
structure of the wind does not change significantly. However, 
the velocity at the base of the outflow is affected. The change 
in velocity due to radiative acceleration is shown in Figureflll 
This figure shows both the poloidal velocity as a function of 
distance along the flowline, and the variation in that velocity 
with iterations of the model, therefore showing the conver- 
gence in the model. Figure ^2 shows that line-driving near 
the disk does significantly accelerate the wind, but magnetic 
driving determines the terminal velocity at larger radii. This 
figure also displays how the code converges; the relatively 
slow convergence during the first four iterations is the result of 
the program slowly increasing the radiative acceleration to the 
computed value (increased slowly in order to avoid severely 
overestimating the radiative acceleration in the lines and caus- 
ing sudden deceleration in later iterations). In the later itera- 
tions, the calculation converges to a final velocity profile using 
the full radiative acceleration. This profile shows the affect 
of line-driving near the base of the wind, where the velocity 
increases above that of the initial magnetocentrifugal wind. 
However, as already mentioned, the magnetocentrifugal wind 



(in this model, where L/LEdd = 0.01) still determines the ve- 
locity at large distances. At those distances, the gas is too 
ionized to be appreciably accelerated by line driving. 

It is important to note that considering non-Sobolev effects 
leads to large changes in Mij nes , the line force multiplier, found 
in the above calculations: "capping" the Sobolev absorption 
length-scale by the actual absorbing column length leads to 
smaller optical depths and larger line accelerations. This is 
illustrated in Figure ^] where the fiducial model has been 
calculated with the Sobolev ap proxim ation as well as our non- 
Sobolev treatment (see Section l3.3.2t . The difference in force 
multiplier is due to the strict Sobolev treatment overestimat- 
ing the column for the low-acceleration gas near the base of 
the wind. This relatively straightforward modification is very 
important to correctly estimate the optical depth, as can be 
seen in Figure [T2l 

Overall, in this section, the fiducial model has shown how 
the wind velocity, number density, column density, and radia- 
tive acceleration all interact to determine the final state of a 
magnetocentrifugal wind. These components have not previ- 
ously been self-consistently combined in a magnetocentrifu- 
gal model, and so yield a new look at the state of these winds. 
In addition, the importance of the number densities and col- 
umn densities to the final result are most apparent, and clearly 
merit further investigation, which will be addressed in ^5] 

5. DEPENDENCE OF WIND STRUCTURE ON MODEL PARAMETERS 

Having analyzed the fiducial model in detail, and observed 
how that model's properties change with height, how sensitive 
are the trends in |4]to those fiducial parameters? This is an 
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FIG. 7. — The geometry, dynamics, and photoionization state of an illuminated, radiatively-accelerated magnetocentrifugal wind. Pane (a) shows the height 
of the wind (in units of the launching radius) with respect to distance along the flowline, s. The arrows in pane (a) give the positions of the columns shown in 
Figure 0(1O 23 , 10 225 , 10 22 , 10 21 cm" 2 ; decreasing with increasing s and z). Pane (b) shows the various components of the velocity: poloidal, radial, azimuthal, 
an vertical, respectively, in units of the Keplerian speed at the base of the flow (for the fiducial model, that base Keplerian speed, v^ o, is 6.65 X 10 3 km s" 1 ). Pane 
(c) displays various densities: the total hydrogen and electron number densities as well as the total hydrogen column and the effective Hydrogen column (defined 
as the column of cold, neutral gas of solar abundances that would produce the same obscuration at 1 keV). All densities are plotted in units of the number density 
or total column density at the base of the wind. Pane (d) shows both the electron temperature in the gas and the ionization parameter, U = n 7 /kh where « 7 is the 
number density of hydrogen-ionizing photons. 



important question, and one of the key attributes of this self- 
similar model is that it allows some flexibility in the selection 
of initial parameters. In this section, we test for variations in 
the wind by examining how the wind changes as parameters 
are modified. 

5.1. Variations with Shielding Column 

One of the most difficult issues for radiative driving in 
AGNs is over-ionization of the wind. As shown in &\ as the 
wind accelerates and its density decreases, the magnetocen- 
trifugal outflow can easily become too ionized to be efficiently 
accelerated to escape velocity solely by atomic lines. This is 
the problem of the "shielding gas" that was mentioned in §l\ 
for pure radiative line-driving, some shielding gas is required 
to intercept the X-ray ionizing radiation so that the remaining 
UV resonant line photons can be absorbed by the wind and 
radiatively accelerate it to the escape velocity. 

Some important papers have already been dedicated 
to examining the concept of shielding gas, such as 
IChelouche & Netzerl ll2003bl considering very detailed pho- 
toionization simulat ions of gas shie l ds wi th constant col- 
umn density) and iProea & Kallmanl (120041 where multidi- 



mensional hydrodynamics simulations with approximate ra- 
diative effects are considered). In contrast, the models pre- 
sented here include a shield where the column density varies 
with height in a shielding wind, and where detailed photoion- 
ization simulations can be employed. 

The MHD wind model presented here already launches 
a wind magnetocentrifugally, so it is immune to con- 
cerns of overionization. Is it therefore possible for a 
magnetocentrifugally-driven wind, with its commensurate 
drop in column density with height above the disk, to act as a 
shield, allowing for more efficient radiative acceleration be- 
yond it? We can test this question by simply varying the 
shielding column in the fiducial model and checking the radia- 
tive acceleration seen by a wind launched behind the shield. 

As shown in Figure [D] as the shielding column is in- 
creased from A^e.shieid.o ~ 10 21 to 10 24 cm" 2 , line-driving in 
the wind increases from Tn nes ~ 0.05 to ri; nes ^0.12. (Recall 
that for this model, continuum driving is r conanuum ~ 0.01.) 
This shows that, for large shielding columns (iVH,shieid,o ~ 
10 24 cm -2 ), line driving can be up to an order of magnitude 
more effective than continuum-driving at the base of the wind. 
This increase in acceleration is due to the absorption of the 
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FIG. 8. — Variation in all components of the velocity between the initial, pure 
driven model. The poloidal velocity (shown in pane a) shows the effects of 
approximately a factor of two for this low L/Z-Edd =0.01), as do the poloidal 
azimuthal component of the velocity is relatively unaffected (pane c). 
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magnetocentrifugal wind model and the final magnetocentrifugal & radiatively- 
radiative acceleration (increasing the velocities near the base of the wind by 
velocity's components, the radial and vertical velocities (panes b and d). The 



ionizing radiation by the shield, which allows a lower ion- 
ization state in the wind, and more line-driving due to more 
atomic lines. It is also apparent that, as the shielding is in- 
creased, the resultant lower total flux at the base of the wind 
means that the onset of significant radiative acceleration is 
delayed: this accounts for the offset of maximum radiative 
acceleration from the disk surface as the shielding column is 
increased. 

Further, in Figure the ratio of radiative acceleration to 
magnetic acceleration along a streamline is shown. Since the 
MHD effects are still supplying most of the acceleration (with 
an acceleration roughly equal to and opposite that of gravity), 
the ratio of radiative to magnetic acceleration looks much like 
that in Figure [D] Magnetic effects dominate in these wind 
models, even when large columns of shielding are included. 

We have therefore shown that a magnetocentrifugal outflow 
can act as a shield and increase the efficiency of line-driving 
in the wind. However, it can also be seen that line-driving is 
important in these models only at the base of the wind. This 
arises not only from the drop in the shield's column density 
with height above the disk, but the drop in the wind's density 
as well (and the commensurate rise of the ionization parame- 
ter). 

5.2. Variations with Initial Density 



Owing to the increase of line-driving with decreasing ion- 
ization parameter, higher accelerations would also be ex- 
pected at higher densities. Thus, we investigate the effect of 
changes in the initial density in the wind in FiguresH31andll6l 

Displaying the effect of a range of initial densities, Fig- 
ure[2]shows that line-driving is only effective in these mag- 
netocentrifugal winds at relatively high densities. Since con- 
tinuum driving is approximately constant (and relatively inde- 
pendent of density) at a/g ~ 0.01 for L/Lndd = 0.01, any line- 
driving below that level is insignificant for these winds. Thus, 
for initial densities «o < 10 9 cm" 3 , line driving falls below the 
level of continuum driving and ceases to be important. For the 
highest density tested, «o = 10 11 cm -3 , line driving dominates 
continuum driving for all locations in the wind. (The varia- 
tions in each acceleration curve as a function of s shown on 
this plot are chiefly due to variations in ionization parameter: 
line-driving is high near the disk due to shielding and the rela- 
tively high density, and rises towards the end of the streamline 
due to the dropping flux levels at large distances.) 

Since no clearly has a great impact on the radiative acceler- 
ation, how do such changes affect observables, such as the ve- 
locity? Figure [K)] shows how the variation in radiative accel- 
eration affects the poloidal velocity (v p ) of the outflow. As the 
initial density and radiative acceleration increase, the wind's 
velocity shows substantial variations from the pure magneto- 
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FIG. 9. — Variation in density between the initial, pure magnetocentrifu- 
gal wind model and the final magnetocentrifugal & radiatively-driven model. 
Shown here is the change in the hydrogen density with height, but the elec- 
tron density and columns change similarly: the added radiative acceleration 
yields a drop in density near the base of the disk, as can be seen in all of the 
various density measurements. 
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FIG. 11. — Poloidal velocity for the fiducial wind as a function of dis- 
tance along the flowlines, showing the evolution of the velocity over several 
iterations of the code. Line driving does significantly affect the magneto- 
centrifugal wind near the base, but magnetocentrifugal driving dominates at 
larger distances for this Eddington ratio. (In the first iteration in this plot, 
only magnetocentrifugal driving has been considered, so the "Iteration #1" 
curve shows the pure magnetocentrifugal wind case.) 
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FIG. 10. — Line radiative acceleration, continuum radiative acceleration, 
and magnetocentrifugal acceleration along the gas streamlines, compared in 
the fiducial model. Note that at the high densities and columns near the base 
of the wind, line-driving dominates continuum driving, but that both are less 
than the magnetic driving for a system at Z./iEdd = 0.01. At larger distances, 
continuum driving dominates line driving in the more highly ionized gas. As 
expected, continuum-driving is approximately of order the Eddington ratio 
for the highly-ionized portion of the outflow. 
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FIG. 12. — Comparison of the line force multiplier, Mij nes (see Eq.l4l. 
for the case of our non-Sobolev method (defined in Section 13.3.21 and the 
strict Sobolev calculation. The over-estimate of the absorbing column in the 
Sobolev approximation leads to higher opacity in the gas (higher r) and thus 
lower acceleration. 



centrifugal model (which dominates the «o = 10 7 cm" 3 and 
«o = 10 s cm" 3 models). In the case of «o = 10 11 cm" 3 , where 
the greatest difference in v p is seen, the velocity increases 
by a factor of ~ 2 to 3 close to the disk. Beyond the re- 
gion close to the disk (i/rrj > 1), however, magnetocentrifu- 
gal driving still dominates the final velocities for these winds. 
But the velocity differences near the disk may be observa- 
tionally important, especially if acceleration near the base of 
the wind is the source o f single-peaked emission lines (as in 
iMurrav & Chiangll997li . If true, such emission lines may be 
critical in testing the differences between wind models. 

5.3. Variations with Radiative Column 



Having already tested the more obvious parameters of the 
initial density and shield column density, we now turn to 
one of the most crucial parameters for the efficiency of line- 
driving: the optical depth in the lines. 

Under normal circumstances, where the gas velocity is suf- 
ficiently low, or where the gas column is very low, the optical 
depths in the lines are simply governed by the ionic columns 
themselves. For large accelerations or large columns, how- 
ever, the optical depths are dominated by the Sobolev length, 
which is defined as the distance over which the relative ve- 
locity between atoms is equal to the thermal width, so that 
a photon emitted by one atom c ould be absorbed by another 
within that Sobolev length (see 33.3. 11 1. With the Eddington 
ratio in this magnetocentrifugal wind model (L/Ledd = 0.01), 
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FIG. 13. — Variations in line-driving as the magnetocentrifugally-launched 
shielding column is increased. With greater shielding columns, the wind is 
less highly ionized, leading to higher radiative acceleration near the base of 
the wind. Also, with increased shielding comes lower flux levels at the base 
of the wind, which displaces the onset of line-driving to larger distances from 
the launch point. (Continuum-driving is approximately constant at a/g ~ 
0.01.) 
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FIG. 14. — As in Figure fT3l but comparing the strength of the total radiative 
acceleration (both line and continuum acceleration) to the acceleration due to 
magnetic fields, where both are normalized relative to gravity. To compare 
with the magnetic field ac celer ation, the radiative acceleration is scaled by a 
geometrical factor (see Eg . 1131 that yields the radiative acceleration along the 
streamline. 



both regimes can be important. Depending on the initial pa- 
rameters prescribed, the wind column can be small enough 
such that the column alone determines the opacity in the wind 
(instead of the velocity) and therefore the amount of observed 
acceleration, so that the Sobolev length is not important. This 
is demonstrated in Figure [n] where the variation in line- 
driving with the radiatively accelerated wind column is pre- 
sented (the shielding column is held constant). For the larger 
columns, the optical depth in the lines increases and the ra- 
diative acceleration decreases. For the smallest columns, very 
large radiative acceleration is predicted due to the low opac- 
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FIG. 15. — Variations in line-driving with changes in the density at the 
base of the wind. Higher densities lead to smaller ionization parameters and 
much larger acceleration. (Continuum-driving is approximately constant at 
a/g ~ 0.01.) 
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FIG. 16. — As in Fig. 1151 except that the change in velocity along the 
streamline is shown for the variety of densities. For this Eddington ratio, 
the magnetocentrifugal wind still dominates, but as the density increases, the 
importance radiative driving increases. 

ity in the lines. This is critical for these models, for at sig- 
nificantly high column, the wind would see no significant line 
driving (this is true for magnetocentrifugal wind columns with 

Akrad.o > 10 22 cm" 2 ). 

5.4. Modifying the Eddington Ratio 

Of key importance to applications to AGN is understanding 
the acceleration of outflows as a function of the Eddington 
ratio, L/LEdd- To investigate the impact of varying Eddington 
ratios in our model, we present Figure^] which displays the 
radiative acceleration (both the combined line and continuum 
acceleration in panel a as well as the line driving in panel 
b) in the fiducial model for three different Eddington ratios: 
L/L^dd = 0.001,0.01, and 0.1. 

The largest variation in FigurelT8h is the continuum driving 
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FIG. 17. — Variations in line-driving as the radiative wind column is in- 
creased. As the wind thickness increases, the optical depth within of the 
lines within the magnetocentrifugal wind increases and line-driving drops in 
strength. (Continuum-driving is approximately constant at a/g ~ 0.01.) 



increasing linearly with the Eddington ratio. As expected, the 
continuum acceleration, relative to gravity, is roughly equal to 
the Eddington ratio. The line driving, on the other hand, can 
be seen in both the deviations from the approximately con- 
stant continuum acceleration in Figure[T8h and in Figure[T8b. 

We begin examining this figure by concentrating on the first 
three models in Figure^] which have -/Vn.shieid.o = 10 23 cm -2 . 
In these models in FigurefT8h. the increase in acceleration due 
to line-driving, relative to the continuum-driving, decreases 
as the Eddington ratio increases. This can also be seen in 
Figure[T8b. where the line-driving peaks near the disk (s / ro < 
0.01) for L/Lem = 0.01, but decreases for Eddington ratios 
an order of magnitude larger (where the gas is overionized) 
and an order of magnitude smaller (where the radiation field 
doesn't have the momentum to accelerate the wind as strongly 
as at L/L Edd = 0.01). 

Now we turn to the fourth model in Figure ED where 
we keep L/Lem = 0.1 but increase the shielding level to 
A^H.shieid.o = 10 24 cm -2 . This model shows the importance of 
shielding gas for L/Lem = 0.1. The increase in L/Lgda and 
the increase in shielding relative to the fiducial model allow 
for increased line radiative acceleration that jumps almost two 
orders of magnitude in strength near the base of the wind. 

6. RESULTS 

If magnetocentrifugal winds power outflows in AGNs, they 
must certainly be affected by the intense radiation field they 
experience; in turn, such winds will also influence the ef- 
ficiency of radiative acceleration. This paper has explored 
the radiative transfer through these magnetocentrifugal winds, 
and how they both are affected by and affect radiative accel- 
eration of outflows from AGNs. The model has been used 
to explore the detailed dynamics and ionization of a fiducial 
magnetocentrifugal disk wind, showing the inter-relation be- 
tween shielding column, initial number density, outflow ve- 
locity, Eddington ratio, and acceleration. 

As a result of this study, these models have shown: 

1. A magnetocentrifugal outflow, acting as a "shield", 
can improve the efficiency of line-driving by factors 



o f app roximately two to three A^j.shieid.o = 10 cm 
( 35. 1> and by up to almo st tw o orders of magnitude 
for A^H.sMeid.o = 10 24 cm" 2 ( 35.4> . A magnetocentrifugal 
wind has the advantage that it can be accelerated with- 
out regard to the ionization state, whereas radiatively- 
driven winds must have a low ionization parameter in 
order for a critical abundance of atomic lines to be 
present. Therefore, magnetocentrifugal winds could 
play an important role in acting as a radiation shield 
and allow large radiative accelerations. It may also be 
possible t hat pressure differences (MCGV95) or disk 
photons (Proga & Kallman 2004) may help "lift" the 
shield; neither of those effects are considered here. 
Later work with this model will include the effect of 
disk-emitted photons. 

2. The efficiency of line-driving is str ongly dependent on 
the density at the base of the wind ( SI5.2> . This is due to 
the very critical dependence of line acceleration on the 
ionization parameter. The lower the ionization state, 
the more lines exist to aid in the momentum transfer 
from outward-streaming photons. The density at the 
base of the disk is therefore crucial to setting to line- 
driving within these magnetocentrifugal models. 

3. Small columns (A^H.rad.o 10 21 cm -2 ) within magneto- 
centrifugal wi nds can be significantly accelerated by 
line-driving ( SI5.3> . This point demonstrates the im- 
portance of "non-Sobolev" effects; i.e., that at low 
columns, the optical depths in the lines drop below the 
opacity given by the Sobolev length, and at such low 
columns, the radiative acceleration can be underesti- 
mated by the simple Sobolev approximation. 

In addition, by examining the fiducial model and the above 
cases where model parameters are varied, these solutions have 
displayed the importance of considering the detailed interac- 
tion between the dynamics and photoionization in AGN out- 
flows. Calculations of the ionization parameter along the flow 
and in the variation of shielding and optical depth along the 
flow are central issues to modeling these winds. The issues 
outlined above are a few of the dependences that arise from 
such modeling, and may indeed (by the variation in accelera- 
tion and therefore velocity along the streamlines) lead to tests 
to observationally determine the physics of wind launching in 
AGNs. 

Thus, while this model has been developed to help address 
the above questions about shielding and the affects of radia- 
tive driving on magnetocentrifugal winds, the solutions avail- 
able are not limited to the above, examined cases. Future 
papers will study further the variation of radiative accelera- 
tion on model parameters such as the SED, atomic line lists 
used to calculate the acceleration, initial densities, Edding- 
ton ratios, and other parameters of the model. The model can 
also be used to explore the absorption and emission features 
from such a wind, as well as, for ins tance, the possible role 
of "clouds" within a continuous wind (Evere ttet alJ2002l) . In 
addition, this model is in no way constrained to only study 
AGNs. The same basic physical framework could also be 
employed to study winds from accretion disks surrounding 
young stellar objects or cataclysmic variables. 
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FIG. 18. — Variation in both continuum and line driving (a) and line driving alone (b) when the Eddington ratio in the fiducial model is modified. The first 
three models all have NH.shield.o = 10 23 cm -2 , while the last model has Wh. shield, = 10 24 cirT 2 . For lower Eddington ratios (L/Lem = 0.001,0.01) the gas is in 
a relatively low ionization state for line-driving to be important, and the acceleration departs from the pure continuum acceleration (which sets the acceleration 
relative to gravity to about the Eddington ratio). For a high Eddington ratio (L/Z-Edd = 0-1). the g as is overionized and only continuum acceleration is important; 
raising the shielding to a nearly Compton-thick iVH.shidd.o = 10 24 cirT 2 lowers the ionization state of the gas and allows much higher radiative acceleration. 
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APPENDIX 

APPENDIX A: DERIVATION OF THE SELF-SIMILAR CENTRIFUGAL WIND EQUATIONS 

In this appendix, a rederivation of the system of self-similar wind equations for the magnetocentrifugal wind is presented. The 
equations utilized in this calculation advance upon those presented in BP82 and KK94: the wind not only has an arbitrary density 
power-law index, b, as in KK94, but energy conservation is not required. Since the radiation field continually inputs energy into 
the outflow, this is an important modification that was not fully considered in the derivation presented in KK94. 

First, a stationary, axisymmetric, ideal, cold MHD flow in cylindrical coordinates (r, <j),z) is assumed. The equations are based 
on both the radial and vertical momentum equations: 



dv r dv r 
Or oz 



v 



d<f> 
r or 



47T 



dB z 
dr 



dB r 
~dz~ 



<9(rfi ) 



4irr dr 



d$ 1 dB 2 1 
P(v-VK=-p---- + -(B.V)* z , 



(Al) 



(A2) 



where v is the fluid velocity and B is the magnetic field. The thermal term is neglected in the limit that thermal affects are much 
less important than magnetocentrifugal and radiative-driving effects. $ is the effective gravitational potential, defined as 

GM. 



$= -[1-1X0)]- 



(A3) 



V + z 2 ) 1/2 ' 

where M, is the mass of the central black hole, and T(9) gives the local radiative-to-gravitational radial acceleration (see eq. (Q). 
These equations are solved by first stipulating mass conservation, 

V-(pv) = 0, (A4) 

and then relating the flow velocity to the magnetic field via 



v(r) = - — — + uj{r) x r, 



(A5) 



47rp(r) 

(e.g., IChan drasekha rll 9~56t iMestelll 1 96 ll) . where k/4ir is the ratio of mass flux to magnetic flux, and uj(r) and p(r) are the field 
angular velocity and gas mass density of the flow, respectively. Both uo and k are constant along magnetic fieldlines. In addition, 
while the specific energy is not constant, the total specific angular momentum 



. rB^ 
l = rv,- — 



(A6) 



is conserved. 
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Self-similarity is then imposed on this system by specifying 

r = [r ax),(/>,roXl (A7) 

v=lC(x)f(x),g(x)J(x)]nfi, (A8) 

where Vk,o is the Keplerian speed at the base of the outflow, v^.o = (GM./ro) 1 / 2 , and the prime indicates differentiation with 
respect to \- At the same time, the above constants are re-expressed in dimensionless form: 

x ^w6^> (A9) 



£a+d) i/2 



vjto, (A10) 



Bpfl 

where B p x> is the poloidal magnetic field strength at the base of the wind. 
As in KK94, a general power-law scaling of the density and magnetic field along the disk's surface is defined: 

pocxro", (All) 

fiooc/ +1)/2 . (A12) 
With this self-similar specification, the radial and vertical momentum equations become, after some simplification: 



k£J £/ s J £ 3 (m-l) 2 SL WJ </ 2 V 2 



(A13) 



(A14) 



where 



£ JS 2 J £ (m-1) 

(A-e 2 ) (-6+1) 
(m-1) 2 

X V(m-l) (m-1) 2 

— (m-/<y+/<xC )=-[i-r(0)]xS +-TT2 ~ + 

Kt,J K^J 1 \ 2 

(x+ ((')(' ("(\~H*> 

w(A-af (t+ "' A : ?2) -f + ^ 

V 2£(m-l) 2 

(m-1) 3 ' 
47T/OV 2 

m = = = square of poloidal Alfven Mach number, (A15) 

p 

&(1 £' 2 )^v 

k = — = dimensionless ratio of mass flux to magnetic flux, (A16) 

B P ,o 

A = r = normalized angular momentum, (A17) 

(GMr )2 

J = i-X^ (A18) 

S=l/Ve + X 2 - (A19) 

(A20) 

The two equations ( IA13i and (IA14> define the differential equations for m' and which are, respectively, the spatial gradient 
in the poloidal Alfven mach number (gradient with respect to height, x) an d the (cylindrical) radial velocity gradient (again with 
respect to x). 

One can see from close inspection of the above equations that many of the terms have a denominator of (m-1), showing that 
when the gas crosses the Alfven point (where m =1), the equations become singular. Rewriting and solving the m' equation for 
the value of m' at the Alfven singular point: 

m' A = 2^[-8xK 2 Am'^7 3 + 4( 1 + b)n 2 X£ 2 J 2 ( X + f f ') + m' 2 ( X + &')(-2k 2 X 2 S+ 

(l+b) + 2k 2 X 3 - 4 X n 2 A? ( A - S)£' + ((1 +b) + 2 x 2 k 2 X(X - SM' 2 + 2k 2 XT(8)J 2 )] / 



(A21) 
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This constraint is used to start the integral at the Alfven point with the value of m' A given by equation ( IA21I . 

As covered in the main text, these equations are solved using a "shooting algorithm," integrating from both the Alfven point 
and the disk surface towards an intermediate point. Matchin g the inte grals of three first-order equations (given by the first-order 
equation for m' and the second-order equation for £" in eas. lA13l and lA14t at the common point allows us to solve for the three 
free parameters in the system: £q, £' a , and xa- 
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